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ABSTRACT 

Selsraologlcal  research  supported  by  Air  Force  contract  F19D28-71-C-0245 
at  the  Lamont “Doherty  Geological  Observatory  Is  summarized  for  the  period 
1  July  1971  to  31  December  1971* 

During  this  period  extensive  analysis  was  done  on  the  systematic  re¬ 
gional  variations  In  stresses.  High  stress  earthquakes  that  usually  ex¬ 
hibit  anomalously  large  hlgh-fhequency  content  are  the  most  likely  ones  to 
be  confused  with  explosion  signals.  Higher  stresses  are  expected  to  be  ac¬ 
cumulated  In  regions  where  new  faults  are  being  created  In  corrpetent  rocks. 

A  possible  temporal  change  of  travel  time  Is  observed  at  some  of  the 
stations  In  western  United  States  that  may  have  been,  at  least  In  part, 
attributed  to  a  change  of  the  stress  conditions  along  the  path. 

The  use  of  40-second,  rather  than  20-second,  surface  wave  anplltudes 
for  nu  -  M  graphs  was  shown  to  have  some  distinct  advantages  for  purposes 
of  detect iSn  ani  discrimination.  The  pronounced  difference  3n  spectral  con¬ 
tent  of  P  waves  from  earthquakes  and  underground  nuclear  explosions  sug¬ 
gests  a  method  for  Identification  that  Is  applicable  for  explosions  large 
eno'i^  to  be  recorded  by  long-period  seismographs  at  teleselsmlc  distances. 
This  method  may  be  particularly  Important  when  long-period  Instruments  are 
disturbed  by  other  events  and  meth^s  of  Identification  that  rely  on 
analysis  of. surface  waves  become  Inapplicable. 


MAJOR  SCIHTOJTC  ACCOMPLISHiyENTS 

In  the  followjng  paragraphs,  scientific  accompllshemnts  are  sumnarlzed 
following  the  itemization  in  the  statement  of  work  of  this  contract. 

Line  -  Item  la 

Operation  of  the  world-wide  Lamont-Doherty  network  of  long-  and  Inter¬ 
mediate-period  seismic  stations  at  the  Palisades  (PAL),  Sterling  Forest 
(SFO),  ani  Ogiensburg  (OGD)  continued  during  the  present  contract  period. 

line  -  Item  lb 


1.  During  the  last  period,  a  study  of  earthquakes  with  anomalously 
large  high-frequency  content  was  begun.  These  events  are  of  special  inter¬ 
est  because  they  are  the  most  likely  ones  to  be  confused  with  explosion 
signals.  High-frequency  earthquakes  are  also  called  high  stress  earth¬ 
quakes  because  high  stresses  in  the  source  region  are  the  most  likely  cause 
for  relatively  strong  high  frequency  content. 

Before  the  discrimination  problem  for  Individual  high  stress  events  Is 
attacked,  w'e  must  recognize  these  events  and  we  may  ask  whether  there  are 
systenw.tlc  regional  variations  In  the  stresses.  Higher  stresses  are  ex¬ 
pected  to  be  accumulated  In  regions  where  new  faults  are  being  created  In 
competent  rocks.  Low  stresses  may  p:’evall  along  well  developed  active 
faults.  The  simplest  method  to  recognize  hl^  stress  events  Is  to  classify 
earthquakes  by  the  ratio  of  a  high  frequency  to  a  low  frequency  amplitude. 

A  more  sophisticated  method  Is  to  obtain  the  body  wave  spectrum  of  an  earth¬ 
quake  over  approximately  two  decades  of  frequency  and  determine  the  stress- 
drop. 


2.  A  paper  by  Moliiar  and  ^ss  in  which  stress-drops  of  shallow  earth¬ 
quakes  in  the  Tonga-Kennadec  arc  were  determined  from  displacement  spectra 
of  body  waves  was  submitted  to  Physics  of  the  Solid  Earth.  They  analyzed 
37  events  and  found  that  earthquakes  with  the  highest  stress-drops  occurred 
outside  the  main  thrusting  zone.  They  concluded  that  the  high  stress-<lrop 
events  (high-frequency  risk)  reflected  deformation  within  one  plate  of 
lithosphere  and  that  these  earthquakes  may  have  created  new  faults.  Molnar 
and  Wyss  studied  18  earthquakes  with  depths  greater  than  normal  In  the 
same  area  and  found  that  earthquakes  of  intermediate  depth  had  higher  stress- 
drops  tnan  sl^low  ones.  This  observation  corroborated  the  results  by 
^ujiura  (1969)  and  Wyss  (1970)  who  both  found  marked  increase  In  higji- 
frequency  content  for  earthquakes  with  Increasing  depth. 

3«  A  paper  entitled  "The  use  of  body-wave  spectra  In  the  determination 
of  seismic  source  parameters"  by  Hanks  and  Wyss  was  submitted  to  the  Bulle¬ 
tin  of  Seismologlcal  Society  of  America.  Teleselsmlc  determination  of  body- 
(P,  S)  spectra,  Interpreted  In  terms  of  the  Brune  (1970)  seismic  source 
model,  are  used  to  estimate  the  parameters,  seismic  moment,  and  source  di¬ 
mension  for  three  large,  shallow,  strike— slip  earthquakes  occurring  on 
nearly  vertical  fault  planes  and  for  which  the  same  parameters  can  be 
determined  from  field  data.  These  earthquakes  are  the  Borrego  Mountain 
California,  earthquake  (April  9,  1968);  the  Mudumu  Valley,  Tbrkey,  earth- 
^2,  1967);  and  the  Dasht-e-Bayaz ,  Iran,  earthquake  (August  31, 
1900;.  A  minimum  estimate  for  the  radiated  energy  Is  low  by  a  factor  of  3- 
10  with  respect  to  the  estimate  obtained  from  energy  -  magnitude  relations 
for  these  three  earthquakes.  The  stress  drops  of  these  events  are  of  the 
order  of  10  bars. 

The  above  studies  underline  the  importance  of  exact  depth  determln- 
atlons  for  more  or  less  shallow  earthquakes.  The  deeper  an  event  the  more 
high  frequencies  will  be  radiated  and  the  more  difficult  will  It  be  to 
discrlmlmte  it  from  explosions.  The  second  conclusion  from  the  above 
studies  ^s  that  in  an  average  sense  there  are  regional  differences  In  the 
high  ft?equency  content.  However,  hlgh-stress  regions  produce  occasionally 
low-stress  earthquakes  and  vice  versa. 

A  paper  entitled  "Seismic  body  waves  In  the  vicinity  of  Mount 
Katmal,  Alaska,  and  evidence  for  the  existence  of  molten  chambers"  by 
Tbs^tu  Matumoto  was  published  In  the  Bulletin  of  the  Geological  Society 
of  America.  The  study  revealed  the  disappearance  of  shear  waves  along 
certain  paths  crossing  the  volcanic  range.  A  close  geographical  correla¬ 
tion  of  these  paths  with  active  volcanoes  strongly  suggests  that  magna  or 
zones  of  partial  melting  exist  and  are  responsible  for  S-wave  shadowing. 

By  the  use  of  the  high  attenuation  of  high-frequency  P  waves,  an  estimate 
of  viscosity  In  a  raagna  chamber  of  10°  cgs  units  is  made. 

5.  A  detailed  mapping  of  the  attenuation  properties  of  the  wedges  of 
mantle  behind  Island  arcs  and  above  the  descending  lithospheric  slabs  Is 
currently  under  way.  In  this  study  pP  and  P  waves  from  Intermediate  and 
deep  earthquakes  as  recorded  on  the  VJWSSN  seismograms  were  used.  Pre¬ 
liminary  results  strongly  confirm  the  existence  of  the  hlgh-attenuation 
zone  behlM  the  Tonga-Kermadec  arc  as  found  by  Barazangi  and  Isacks  (1971). 

A  study  of  all  the  Earth’s  Island  arcs  and  Islaiid  arc-llke  structures  Is 
being  can’led  out. 


3. 


6.  An  algorithm  has  been  developed  to  calculate  difference  formulae 
for  boundary  and  Interface  conditions  arising  In  the  solution  of  the 
elastic  wave  equation  by  finite  difference  methods.  Using  the  calculated 
formula  for  the  comer  of  a  quarter  space,  a  Rayleigh  wave  was  propagated 
around  the  comer.  Agreement  with  model  experiments  was  within  the  com¬ 
putational  error.  A  similar  experiment  was  performed  for  a  three-quarter 
space  and  agreement  was  obtained  with  the  few  published  values.  The 
algorithm  can  be  used  for  any  model  whose  boundaries  and  Interfaces  can 
be  approximated  by  straight  lines  such  as  wedges,  Imbedded  material,  sur¬ 
face  elevations,  etc. 

7.  The  elastic  waves  radiated  by  a  stress-free  surface,  which  grows 
within  a  pre-stressed  medium  was  made.  Earthquakes  may  reasonably  be 
modeled  In  this  way,  and,  for  cases  with  an  elliptic  fault  plane,  the 
radiated  waves  can  (surprisingly)  be  calculated  by  Cagnlard-de  Hoop 
theory.  A  simple  method  has  also  been  found  to  obtain  the  far-fleld 
pulse  shape  from  given  motions  on  the  fault  plane. 

8.  A  study  Is  being  conducted  to  find  the  Initial  phase  patterns  as¬ 
sociated  with  azimuth  for  nuclear  explosions  at  NTS  and  nearby  natural 
events.  Long-period  WWSSN  Raylel^  wave  data  Is  used  to  calculate  the 
initial  phase  patterns.  It  Is  hoped  that  a  specific  phase  relation  can  be 
used  In  discrimination-detection  work. 

Line  -  Item  Ic 


1.  A  paper  of  focal  mechanisms  for  shallow  earthquakes  occurring  from 
196^  to  June,  1969,  In  the  region  fron  New  Guinea  to  the  New  Hebrides 
titled  "Focal  mechanisms  and  plate  tectonics  of  the  southwest  Pacific”  by 
Johnson  and  Molnar  was  submitted  to  the  Journal  of  Geophysical  Research. 
Ninety-six  new  focal  mechanisms  were  detennlned  for  earthquakes  on  the 
belt  of  seismic  activity  separating  the  Pacific  and  Australian  plates.  The 
direction  of  convergence  of  these  plates  varies  from  NE-SW  to  east-west. 

The  Australian  plate  underthrusts  the  Pacific  plate  to  the  ENE  under  the 
Solomon  and  New  Hebrides  Islands,  and  overthrusts  the  Pacific  fr  the  east 
along  the  Tonga-Kermadec  arc  and  North  Island  of  New  Zealand.  Between  the 
southern  part  of  the  New  Hebrides  arc  and  the  northern  end  of  the  Tonga 
arc,  several  mechanisms  are  consistent  with  a  zone  of  transform  faulting 
between  the  arcs,  but  the  tectonic  boundary  In  this  area  Is  clearly  more 
caiplex  than  a  simple  transform  fault.  The  data  for  the  Macquarie  Ridge 
concur  with  the  Idea  that  the  pole  of  rotation  for  the  Pacific  and  Australian 
plates  Is  nearby  and  to  the  east  of  this  feature.  The  data  also  suggest  a 
NNE-SSW  convergence  of  the  Pacific  plate  and  the  Australian  plate  In  north¬ 
western  New  Guinea. 

2.  A  paper  entitled  "Earthquake  fault  parameters  and  tectonics  In 
Africa"  by  ^!^sha  and  Molnar  has  been  submitted  to  the  Journal  of  Geophysi¬ 
cal  Research.  Fault  plane  solutions  of  earthquakes  In  southern  Africa  Indi¬ 
cate  that  the  least  compressive  stress  is  oriented  approximately  east- 

west  nearly  parallel  to  that  in  the  northern  part  of  the  rift  system.  Seis¬ 
mic  moments,  source  dimensions,  and  stress  drops  were  determined  fur  el^t 
earthquakes  frcan  body  and  surface  wave  spectra  using  the  theory  of  Brune 
(1970).  Spectral  estimates  of  these  quantities  for  the  I966  earthquake  In 
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the  Republic  of  Zaire  agree  well  with  those  observed  in  the  field.  Relatively 
higher  stress  drops  are  found  for  events  not  associated  with  rift  faulting. 

If  higher  stress  drops  indicate  higher  rock  strength,  then  these, 'with  other  ' 
geological  and  geophysical  data,  suggest  tliat  the  northern  part  of  the  rift 
system  is  similar  to  ocean  ridges  and  behaves  as  a  plate  boundary,:  but  the 
southern  part  is  different  and  is  not  a  plate  bouhdfery.  The.  tectonics 
associated  with  the  northern  part  of  the  rift  system  appear  to  migrate  ■ 
southwards .  ' 

,  I 
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A  paper  entitled  "Mantle  wave  analysis  by  a  phase-equallzatlon-pid-?um- 
method  for  the  ftontana  LASA  long-period  data"  by  Kazuo  Hamada  was  published 
in  the  Bulletin  of  the  Selsmologlcal  Society  of  Amerlba.  A  paper  entitled 
"Mantle  Rayleigh  waves  for  shield,  oceanic,  and  tectonic  areas  uslpg  LASA 
long-peri.od  data"  by  Kazuo  Hamada  was  submitted  to  the  Journal  of  Geophysical 
Research. 

By  using  LASA  long-period  digital  data  and  high-galh  digital  data  at 
Ogdens burg.  New  Jersey,  group  velocities  of  mantle  Love  and  Rayleigh  waves 
were  determined  by  mieans  of  a  recently  developed  new  methods  Combining  these 
group  velocities  with  Kanamorl's  phase  velocity  data,  regionalized  shear- 
velocity  models  of  the  upper  mantle  were  inferred.  Important  conclusions 
are:  a.  Shear  velocity  differences  among  the  different  tectonic  regions 
exist  at  depths  less  than  300  km.  b.  Shear  velocities  In  the  upper  mantle 
are  obviously  higher  for  the  shield  areas  and  lower  for  the  tectonic  areas 
than  those  for  the  oceanic  areas;  the  velocity  difference  between  the  shield  , 
and  tectonic  models  Is  approximately  10%  at  depths  less  tLan  100  km,  de¬ 
creasing  gradual.ly;  5%  at  a  depth  of  200  km.  c.  Ih^  shield  data  do  not 
require  a  strong  low  velocity  channel  In  the  upper  mantle  like  the  oceanic 
model,  d.  Ihe  tectonic  data  do  not  necessarily ‘require  the  presence  of 
high  velocity  lid  Just  beneath  M-dlscontlnulty .  e.  The  oceanic  data  re¬ 
quire  obviously  a  strong  low  velocity  channel  In  contrast  vd.th  the  other  tvro ' 
models,  f.  The  dispersion.. of  Love  and  Rayleigh  v/aves  for  the  tectonic 
regions  can  be  explained  well  by  the  different  shear  Velocity  structures  ' 
rather  than  by  the  same  one;  the  shear  velocities  expected  from  Inve  waves 
are  higher  than  those  expected  from  Rayleigh  waves  by  0.2  km/sec.  The 
velocity  differences  are  concentrated  at  depths  from  ,150  to  300  kmi.  This 
discrepancy  requires  anisotropy  or  an  equivalent  laminar-melting  structure.  ' 

I 
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I 

1.  A  paper  entitled  "Microseismicity  and  tectonics  of .the  Nevada  seis¬ 
mic  zone"  by  Prank  J.  Gumper  and  Christopher  Scholz  was  submitted  to  the  , 
Bulletin  of  the  Selsmologlcal  Society  of  America.  Microselsmlclty,  com¬ 
posite  focal-mechanism  solutions,  and  previously  published  focal  parameter 
data  are  used  to  determine  the  current  tectonic  activity  of  the  prominent 
zone  of  seismicity  In  western  Nevada  and  eastern  California,  termed  the 
Nevada  Seismic  Zone.  The  microselsmlclty  substantially  agrees  ^-ri-th  the 
historic  seismicity  and cfellneates  a  narrow,  major  zone  of  activity  that  ex¬ 
tends  from  Owens  Valley,  California,  north  past  Dixie  Valley,  Nevada.  Focal 
parameters  indicate  that  a  regional  pattern  of  NW-SE  tqnsldn  exists  for  the 
western  Basin  and  Range  and  Is  now  producing  crustail  extension  within  the 
Nevada.  Seismic  Zone.  An  eastward  shift  of  the  sqlsmlc  zone  along  the. Excel¬ 
sior  Mountains  and  left-lateral  strike-slip  faulting  determined  from  a  com- 
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poslte  focal  mechanism  Indicate  transform  type  faulting  between  Mono 
lake  and  Pilot  Mountain.  Based  on  these  results  and  other  data,  it  is  sug¬ 
gested  that  the  Nevada  Seismic  ^ne  is  caused  by  the  interaction  of  a  west¬ 
ward  flow  of  mantle  mater^lal  beneath  the  Basin'  and  Range  Province  with  the 
boundary  of  the  Sierra ,  Nevada  bathollth. 

I  2.  A  paper  entitled  VP-wave  spectra  fr'-om  underground  nuclear  explo¬ 
sions”  by  Peter  Molnar  was  submitted  to  the  Geophysical  Journal  of  Royal 
■  Astronomical  Society,  Inndon.  This  paper  reports  observations  of  P-wave 
spectra  from  underground  nucleap  explosions  for  periods  greater  than  about 
1,  s.  The  spectra  have  a  maximum  at  about  2-3  ■  s  and  decrease  rapidly  at 
longer  periods.  Two  phenomena  could  cause  such  a  decrease:  (1)  The  sur¬ 
face  reflection  (pP),  with  l|ts.  polarity  opposite  to  that  of  the  direct  P 
wave,  might  (differentiate  the  long-perlcxi  signal  and  cause  a  modulation  of 
.  2|sinajto/2|  vher^  tQ  is  the  Interval  of  time  between  P  and  pP;  or  (2)  the 
soiree  time  function  for  the  pressure  on  the  boundary  of  the  elastic  zone 
suiroundlng  the  explosion  ml^t  more  neai^ly  resemble  an  Impulse  than  a  step 
'  f\jnctlon  as  often  assumed.  The  data  are  consistent  with  the  simultaneous 
occurrence  of  both  phenomena,  bub  they  are  of  insufficient  quality  or 
quantity  to  establish  which  is  more  Ijiportant.  If  other  information  shows 
that  the  reflection  does  not  have  a  larg?  effect  on  the' observed  spectrum, 
then  the  data  show  that  the  source  time  function  is  primarily  an  Impulse . 

'  I  .  ' 

In  addition,  the  pronounced  difference  in  spectral  content  of  P 
waves  from  earthqual<es  and  underground  nuclear  explosions  suggests  a  method 
for  identification  that  is  applicable  for  explosions  large  enough  to  be 
recorded  by  long-period  seismographs  at  teleselsmlc  distances.  This  method 
may  be  particularly  important  when  long-period  Instruments  are  distrubed  by 
'  other  events  and  methods  of  identification  that  rely  on  analysis  of  surface 
waves  become  inapplicable.  , 

I  !  , 

3*  A  paper  entlt]|.ed  "Excitation  of  seismic  surface  waves  with  perlexis 
of '15  tbiTO- secor/is  for  earthquakes  'and  underground  explosions"  by  J.  Savlno, 
L.  R.  Sykes,  C.  Llebermann,  and  P.  Molnar  was  published  in  the  Journal 
of  ^ophyslcal  Research.!  Using  ciata  from  hlgh-galn  seismographs  operating 
^  in  the  deep  (543-meter)  mine  observatory  at  Ogdensburg,  New  Jersey,  we 
analyzed  the  excitation  of  long-perlbd  surfa,ce  waves  by  eai’thquakes  and 
underground  explosions  for  four  different  regions  of  the  world:  western 
Ihited  States,  the  Aleutians,  Novaya  Zemlya,  and  central  Asia.  The  most 
lii?)ortant  result  of  this  study  is  that  dljscrlminatlon  between  earthquakes 
and  undergrpund,  explosions  on  an  Mg-^rib  basis  is  enhanced  when  the  ampli¬ 
tudes  (pf  Rayleigh  waves  with  periods i  near  4o  sec,i  rather  than  20  sec,  are 
used  to  determine  Mg.  A  discriminant  based  on  surface  waves  with  periods 
near  40  sec  is  particularly  advantageous  because  of  a  very  stable  and  pro¬ 
nounced  minimum  in  earth  noise  in  this  period  range.  Love  waves  also 
yield  a  distinct  separation  between  earthquakes  and  explosions  on  the  basis 
of  Mg  (love)  at  20  and  40  sec  with  m^,  and  can  be  used  for  discrimination 
of  small  events  ifor  v;hlch  no  Rayleigh  vjaves  are  recorded.  The  discrimina¬ 
tion  threshold  at  Ogdensbui^g  is  mh"  =  3.8  for  the  western  United  States  at 
a  distance  of  about  30°  and  mu  =  4.4  for'  the  Aleutians  at  70°.  All  the 
events  studied  could  be  discriminated  by  using  40-sec  Rayleigh  waves. 
Anplltude  spectra  of  Rayleigh  waves  in  the  period  range  15  to  70  sec  were 
computed  for  events  in  the,  western  IMted  States, ’the  Novaya  Zemlya  region. 
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and  the  Aleutians.  Tne  spectral  amplitudes  decrease  more  rapidly  vdth  in¬ 
creasing  period  for  explosions  than  for  earthquakes  in  the  western  United 
States  and  Novaya  Zemlya.  These  spectral,  differences  are  consistent  with 
an  impulsive  source  time  function  for  explosions.  In  the  Aleutians,  how¬ 
ever,  the  shape  of  the  spectrum  of  long-period  waves  for  the  underground 
explosion  Milrow  is  similar  to  spectral  shapes  for  some,  but  not  most, 
nearby  earthquakes.  Two  plausible  explanations  for  this  spectral  sim¬ 
ilarity  are  contamination  of  ohe  Milrow  spectrum  by  earthquake-generated 
Raylel^  waves,  or  the  effects  of  focal  depth  and  radiation  pattern  on  the 
earthquake  spectra. 

A  study  on  a  temporal  change  of  the  travel  time  observed  from  the 
underground  nuclear  explosions  has  been  continued  to  utilize  the  accurately 
established  location  and  origin  time  of  these  events.  The  travel  times 
measured  at  a  station  for  the  repeated  explosions  from  the  Nevada  Test 
Site  show  a  fluctuation  of  as  great  as  over  1  second.  This  fluctuation  of 
the  travel  time  is  primarily  resulted  from  an  emergent  onset  of  the  wave 
and  from  unproper  Identification  of  phases.  By  eliminating  some  of  the 
events  with  smaller  yield,  the  fluctuation  is  greatly  reduced.  Difference 
of  the  distance,  depth,  aM  subsoil  condition  near  the  station  also  contri¬ 
bute  to  the  fluctuation  of  the  travel  time.  Some  of  the  hl^  quality  sta¬ 
tions  (with  low  ground  noise  and  with  accurate  time  marks)  show  consistent 
travel  time,  ani  their  temporal  changes  can  be  examined.  At  the  stations 
MHC(Mt.  Hamilton,  California)  and  KNOT  (Kanab,  Utah),  it  is  observed  that 
the  travel  time  decreases  as  the  seismic  activity  near  the  station  in¬ 
creases  and  the  travel  time  Increases  rather  abruptly  when  the  seismic 
activity  halted.  At  the  most  of  the  stations,  however,  the  relation 
between  the  localized  seismic  activity  and  travel  time  is  rather  compli¬ 
cated.  This  complication  implies  that  the  physical  condition  at  inter¬ 
mediate  depths  that  has  not  been  reflected  to  the  seismic  activity  at 
shallower  depths  nay  play  an  important  role  on  the  temporal  change  of  the 
travel  time. 

Line  -  Item  If 


1.  A  paper  entitled  "Three-dimensional  seismic  ray  tracing  in  a  later¬ 
ally  heterogeneous  spherical  earth"  by  KLaus  H.  Jacob  was  published  in  the 
Journal  of  Geophysical  Research.  Recent  seismologlcal  studies  suggest 
lateral  Inhomogeneities  in  P  and  S  velocities  of  the  mantle  that  are  as¬ 
sociated  with  slabs  of  mobile  lithosphere  descending  into  the  mantle  be¬ 
neath  Island  arcs.  In  special  cases,  travel  times  of  P  traversing  such 
zones  can  differ  by  as  much  as  5  sec  and  of  S  by  up  to  10  sec  from  standard 
travel  times.  In  addition,  such  zones  are  characterized  by  relatively  low 
attenuation  of  S-wave  energy  compared  with  high  attenuation  in  a  broad  zone 
on  the  landward  side  of  the  active  volcanoes.  To  explain  the  observed 
anomalous  travel  times  and  attenuation  phenanena,  it  is  necessary  to  trace 
the  path  of  body  waves  through  laterally  heterogeneous  earth  models .  The 
technique  of  ray  tracing  developed  here  uses  Fennat ’ s  principle  to  obtain 
the  differential  equation  of  a  ray  in  spherical  coordinates.  The  position, 
direction,  and  travel  time  of  the  seismic  v/ave  front  at  any  point  along  the 
curved  ray  path  are  obtained  by  numisrical  integration  of  the  differential 
equation  for  an  assumed  three-dimensional,  continuous  velocity  distribution. 
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Bie  problem  of  representing  a  realistic  three-dimensional  velocity  struc¬ 
ture  in  the  earth  is  solved  in  .a  way  that  is  especially  suitable  for  use 
on  computers.  Some  exarples  for  rays  traversing  an  island-arc  structure 
are  presented.  The  Implications  of  this  method  of  tracing  rays  in  a  lat¬ 
erally  heterogeneous  earth  are  discussed  with  respect  to  seismic  travel- 
time  studies,  interpretation  of  residuals  in  terms  of  tectonic  hetero¬ 
geneities,  source  bias,  and  the  precise  location  of  earthquakes  and 
nuclear  explosions;  dT/dA  measurements  fhom  large  seismic  arrays  and  their 
inversion  to  obtain  details  of  the  velocity  structure  in  the  upper  mantle 
are  also  discussed. 

2.  A  paper  entitled  "Seismic  waves  reflected  from  velocity  gradient 
anomalies  within  the  Earth's  upper  mantle"  by  Paul  G.  PUchards  has  been 
submitted  to  Zeltschrlft  fUr  Geophysik.  Classical  Thomson-Haskell  methods 
have  recently  been  extended,  to  obtain  the  asymptotic  wave  solution  in  a 
stratified  elastic  medium  which  has  both  first  and  second  order  dis¬ 
continuities  in  the  elastic  parameters.  These  methods  are  used  here  in  a 
discussion  of  the  observed  precursors  to  seismic  waves  P'P'.  The  fre¬ 
quency-dependent  reflection  coefficient  R  (=  reflected/incident  displace¬ 
ment  amplitudes)  is  calculated  for  several  models  of  transition  regions 
in  the  Earth's  mantle.  To  generate  observable  precursors  to  P'P',  by 
reflection  from  horizontal  layering  within  the  mantle,  the  thickness  L 
of  the  region  of  transition  is  shown  to  be  much  smaller  than  has  generally 
been  supposed.  This  result  follows  from  the  rapid  decrease  in  R  as  the 
transition  thickness  increases  from  zero  to  one  wavelength.  For  example, 
R(1  second)  >  2  1/2  %  only  if  L  <  4  km.,  even  in  cases  of  1Q%  total 
changes  in  velocity. 
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The  problem  of  representing  a  realistic  three-dimensional  velocity  struc¬ 
ture  in  the  earth  is  solved  in  .-a  way  that  is  especially  suitable  for  use 
on  computers.  Some  examples  for  rays  traversing  an  island-arc  structure 
are  presented.  The  Implications  of  this  method  of  tracing  rays  in  a  lat¬ 
erally  heterogeneous  earth  are  discussed  vdth  respect  to  seismic  travel- 
time  studies,  interpretation  of  residuals  in  terms  of  tectonic  hetero¬ 
geneities,  source  bias,  and  the  precise  location  of  earthquakes  and 
nuclear  explosions;  dT/dA  measurements  from  large  seismic  arrays  and  their 
inversion  to  obtain  details  of  the  velocity  structure  in  the  upper  mantle 
are  also  discussed. 

2.  A  paper  entitled  "Seismic  waves  reflected  from  velocity  gradient 
ananalles  within  the  Earth's  upper  mantle"  by  Paul  G.  Richards  has  been 
submitted  to  Zeltschrlft  fUr  Geophysik.  Classical  Thomson-Haskell  methods 
have  recently  been  extended,  to  obtain  the  asymptotic  wave  solution  in  a 
stratified  elastic  medium  which  has  both  first  and  second  order  dis¬ 
continuities  in  the  elastic  parameters.  These  methods  are  used  here  in  a 
discussion  of  the  observed  precursors  to  seismic  waves  P'P'.  The  fre¬ 
quency-dependent  reflection  coefficient  R  (=  reflected/incldent  displace¬ 
ment  amplitudes)  Is  calculated  for  several  models  of  transition  regions 
Iri  the  Earth's  mantle.  Tb  generate  observable  precursors  to  P'P',  by 
reflection  fran  horizontal  layering  within  the  mantle,  the  thickness  L 
of  the  region  of  transition  Is  shovm  to  be  much  smaller  than  has  generally 
been  supposed.  This  result  follows  from  the  rapid  decrease  in  R  as  the 
transition  thickness  Increases  from  zero  to  one  wavelength.  For  exanple, 
R(1  second)  >  2  1/2  %  only  if  L  <  4  km.,  even  In  cases  of  10%  total 
changes  in  velocity. 
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